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ABSTRACT 

Recent observations and hydrodynamical simulations of star formation inside a giant molecu- 
lar cloud have revealed that, within a star forming region, stars do not form evenly distributed 
throughout this region, but rather in small sub-clumps. It is generally believed that these sub- 
clumps merge and form a young star cluster. The time-scale of this merging process is crucial 
for the evolution and the possible survival of the final star cluster. The key issue is whether 
this merging process happens faster than the time needed to remove the residual gas of the 
cloud. A merging time-scale shorter than the gas-removal time would enhance the survival 
chances of the resulting star cluster. In this paper we show by means of numerical simulations 
that the time-scale of the merging is indeed very fast. Depending on the details of the initial 
sub-clump distribution, the merging may occur before the gas is expelled from the newly- 
formed cluster either via supernovae or the winds from massive stars. Our simulations further 
show that the resulting merger-objects have a higher effective star formation efficiency than 
the overall star forming region and confirm the results that mass-segregated sub-clumps form 
mass-segregated merger-objects. 
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The formation and survival of young, embedded star clusters is an 
important issue not only for the star cluster community but also has 
implications for our understanding of the formation and evolution 
of galaxies as a whole. Stars are now believed to form mostly in 
star clusters, which later dissolve and distribute their stars to their 
host galaxy. 

In recent years, considerable progress in our understand- 
ing of the star formation process has been made and a general 
picture is emerging. Star clu sters form out of coll apsing cloud 
clumps in molecular clouds iTillev & Pudritj [ 20 04). These col- 
lapses a re triggered by turbulent f ragmentation of clouds and their 
clumps (MacLow & Klessen 2004). The clumps are observed to 
be aligned in filamentary structures that can be reproduced by 
supersonic turbulent simulatio n s and they contain many cores 
dBurkert & BodenheimeJ l2000l : iKlessen & Burkertl 1200 if) . Each 
core forms a single star or a binary. The mass function found for 
these cor es is essentially the same as the initial mass function of 
the stars Jjohnstone et"aUl2000l : iGoodwin etaT1l2008h . The cores 
are themselves cl umped in small sub-clumps containing a few to a 
few tens of stars teonnell. Bate & Vindl2003h . While the star for- 
mation efficiency (SFE; i.e. the fraction of gas which ends up in 
the star(s)) in these cores is high, the overall SFE measured over 
the whole molecular cloud is very low, of the order of a few per 



cent dClark & Bonnei]|2004l) . and ^ 40 per cent in clust er-forming 
clumps (i.e. embedded star clusters. [Lada & Lada 2003)). 

The remaining gas does not stay in the newborn star cluster 
but is driven outwards by stellar feedback. In embedded clusters 
containing more than a few hundred stars the feedback consists of 
photo-ionising radiation, the winds of high- mass stars and finally 
the on-set of the first supernova explosions jGoodwirJI 19971) . For 
such clusters the feedback energy can easily be sufficient to un- 
bind the gas leading to a gas expulsion phase whic h is rather short , 
comparable to the crossing time of the star cluster I Rroupa| [2005h . 
Pictures of young massive star clusters, fo r example in the centra l 
region of the Antennae (NGC 4038/4039 Iwhitmore et allll999l) . 
reveal that they are already surrounded by Ha bubbles stemming 
from the gas which was blown out of the star cluster. It has been 
shown that these star clusters can be as young as 5-6 Myr. The out- 
flow velocities of the g as have been measured to be 25-30 km s 
(Whit moreetall l999). which corresponds to gas-evacuation times 
of 0.2 Myr for cluster radii of 4 pc. This is comparable to the cross- 
ing time of a 10 s Mq cluster. As a result of this strong and rapid 
mass-loss the star cluster is left out of virial equilibrium. The ve- 
locities of the stars are too high for the reduced mass of the star 
cluster. Hence, even more mass is lost when stars escape from the 
star cluster. This may finally lead to the co mplete dissolution of 
the cluster (also called 'infant mortality', e.g. lde Griis & Goodwin! 
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Table 1. Table of the initial parameters of our simulations. The first line shows the parameters of our standard model. The columns give the filling factor a, 
the star formation efficiency e and the number of sub-clumps No, followed by the Plummer radius, the cut-off radius, the total mass, the crossing time and the 
velocity dispersion of the star forming region. The next two columns denote the mass in stars and the mass in gas (modelled as analytical background). Finally 
we show the Plummer radius, the cut-off radius, the mass, the crossing time and the velocity dispersion of each single clump. 
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l200ct iGieles & Bastian|[20o3) . But if the star formation efficiency 
is as high as 33 per cent or above, a small bound core remains 
dGoodwid [l997l : (Adamsl [2OO0I : iBoilv & Kroupal [2f3o3al lbh which 
may account for the survival of l ow-mass, Pleiades - type c lusters 
( Kroup a. Aarseth & Hurley 200 3). iGever & Burked [ 2001 ) argue 
that the SFE has to be larger than 50 per cent to get a bound core 
but if the stars have almost no initial velocity dispersion then a SFE 
of 10 per cent could suffice. This means that the measurement of 
the initial velocity distribution in newborn embedded star clusters 
will be crucial to find out whether, and how, star clusters survive. 

All these authors imply that a spherical embedded clus- 
ter has already formed before the remaining gas is blown out 
of the star forming region. In this paper we want to investi- 
gate if this assumption is valid. We start at the stage when the 
sub-clumps have formed (at the end-point of hydrodynamical 
SPH-simulations), but the stars are still in their formation pro- 
cess and therefore have not started to blow away the remaining 
gas. Once the stars start to shine, their stellar winds blow away 
the residual gas. This gas-removal happens on very short time- 
scales, but is longer fo r more massive gas-cores at a given SFE 
jParmentier et alj l2008h . For clusters with masses of less than 
10 5 Mq, this gas-removal tim e is shorter than the crossing-time 
of th e embedded star cluster dBaumgardt. Kroupa & Parmentied 
2008). For massive star clusters, the explosion of the first su- 
pernova will remove the residual gas. This is in principle a dif- 
ferent an s atz th a n most of the numer ic al studies so far (e.g. 
Goodwirj [l998l: IScallv & Clarke! |2002|; iGoodwin & Whitword 
2004L [McMillan. Vesperini & Portegies Zwartl l2007t) . which start 



with a (mostly cold) clumpy structure without the gaseous back- 
ground. 



In this paper we try to answer the following, related, questions: 
Is there sufficient time before the gas is driven out for the sub- 
clumps to merge and form a spherical, embedded cluster? If so, 
what is the effective SFE of these embedded clusters? By means of 
numerical simulations we investigate the time-scale of this merging 
process. 



2 SETUP 

The setup of our simulations is based on a typical outcome of hy- 
drodynamical si mulations of a star forming region inside a giant 
molecular cloud dBonnell. Bate & Vine 2003]). We model the star 
forming region as a Plummer sphere with a mass of M sc = M = 
1000 Mq and a Plummer radius iJ^j = 0.2 pc. We assume that in 
this cloud stars have formed in small sub-clumps with an overall 
star formation efficiency e — 0.32. This gives us 320 Mq in stars, 
which we divide into No = 16 sub-clumps of M c \ — 20 Mq. 
Each of these sub-clumps is modelled as a small Plummer sphere 
with a Plummer radius of R p i — 0.01 pc and a cut-off radius of 
Rcut = 0.05 pc. These clumps are distributed themselves accord- 
ing to the Plummer distribution function of the whole gas cloud 
taking a cut-off radius of 1 pc. Because we do not consider the hy- 
drodynamics of the gas, the remaining mass in gas is modelled as a 
time- varying analytical background Plummer potential. The poten- 
tial is kept constant except when stated otherwise. 

To put the results of this model (our standard model) into con- 
text we also perform a parameter study, in which we vary R^, No 
and e. For clarity we give the parameters of our standard model as 
well as for the comparison runs in Tab.[T] 
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The study of the merging behaviour of these sub-clumps suf- 
fers heavily from the low number statistics (No = 16). For each set 
of parameters, we therefore run up to 5 simulations starting from 
different random realisations and calculate the average of the re- 
sults. 



W e use the particle-mesh code SUPERBOX (Fellhau er et al.l 
l200ti) to perform the simulations. This code has the advantage that 
we can use an arbitrarily high particle number to model each sub- 
clump, because the particles in a collision-less code like SUPER- 
BOX represent phase-space elements rather than single stars. We 
therefore model each clump with 100,000 particles. Furthermore 
SUPERBOX offers two levels of high resolution sub-grids which 
stay focused on each object (i.e. sub-clump) and give us the reso- 
lution needed to follow the merging process in detail and to access 
the time-scales. 

As with all particle-mesh codes, developed to simulate galax- 
ies, the effects of two-body relaxation in our code are almost com- 
pletely suppressed. However, we note that the internal relaxation 
times of our sub-clumps is rather low and of the order of a few 
kyr (7300 yr in the standard model). Although this implies that 
the life-times of those sub-clumps are very short, they will still 
be larger than the time needed to have a first passage through the 
centre of our star forming region, i.e. the time needed to have in- 
teractions with other clumps and the central merger object. These 
interactions are a very violent process which will either destroy the 
clumps and distribute their stars in the central region or reset their 
internal clocks. Since internal relaxation will tend to speed up the 
dissolution of a star cluster, by neglecting internal processes we are 
therefore prolonging the merging process rather than shortening it. 

A study by ISensui. Funato & Makino! j2OO0t) showed that the 
internal structure of galaxies does not play a role in the merging 
time-scales (inside a galaxy cluster) - only the distribution of galax- 
ies inside the cluster matters. This argument should also hold for 
star clusters in a star cluster complex and we believe we can also 
apply it to the problem of this paper. The key quantities are the dis- 
tribution of the relative velocities of the sub-clumps within the star- 
forming region and how much energy they can absorb with each 
interaction. These quantities depend only on the total masses and 
effective radii of the sub-clumps, and hence are reasonably mod- 
elled in our simulations. 

Dissolving sub-clumps (due to the fast two-body relaxation) 
will distribute their stars preferentially in the centre, i.e. in and 
around the merger object. This process would also tend to reduce 
the time-scale for star cluster formation rather than prolonging it. 
As before, our simulations are therefore conservative as they pro- 
vide an upper limit to the total amount of time required to form the 
merger object. 

We are aware that there are also a lot of stars which are ejected 
from the clumps with high velocities due to close encounters within 
the sub-clumps and which are subsequently lost to the forming star 
cluster. However, these stars only constitute a small fraction of the 
total stellar content of the sub-clumps. 

We are therefore confident that the conclusions from our sim- 
plified approach to this problem will also hold for more sophis- 
ticated models, including two-body relaxation, which we plan to 
carry out in the future. 

T o make the comparison with previous work iFellhauer et al.l 
l2002b easier we try to adopt the same notation of dimensionless 
parameters and define as a measure of the concentration of the sub- 



clumps in the gas cloud (filling factor) 

rjclump 



R 



(l) 



pi 



We omit the index clump in the following text. We furthermore de- 
fine a dimensionless time as 



t 

T'sc ' 
cr 

where is the crossing time of the star forming cloud. 



3 THEORETICAL BACKGROUND 
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A theo ry of merging star clusters was developed by Fellhau er et al.l 
d2002h . They investigated the merging time-scales and efficiencies 
of young massive star clusters forming in confined areas which they 
called star cluster complexes. In their case, the star clusters were 
merging in their own potential and not in a background potential 
of the residual gas. Here we want to investigate if this theory also 
applies in the case of stellar sub-clumps forming in the potential of 
the residual gas inside a giant molecular cloud. 

In t his section we recap itulate their theory and adapt it to our 
problem. Fellha uer et al.l J2002I) esta blished first a velocity c riterion 
for merging based on the work of iGerhard & Fail Jl983h which 
reads 

GM pl 



1 2 

2 tn> 
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pi 



(3) 



The right-hand side is a value of the order of the internal veloc- 
ity dispersion of our sub-clumps, while the left-hand side is of the 
order of the clump-velocity dispersion in the star forming region. 
Only if this formula is satisfied should merging between objects be 
possible. If we insert now v typ w \/2o"pi and the definitions for a 
and e we derive 



(4) 



N a 32 
e ^ 37T 

However, Fellha uer et al ] l l2002h noted that in simulations 
with large values of a the right-hand side of equation [3] has to be 
replaced by the unknown values of the merger object. Thus, equa- 
tion [3] gives us a threshold for low-mass and low-a simulations, 
whil e in cases with a high fi lling factor this criterion fails. 

iFellhaueret all (2002) identified two regimes of the merging 
process. As long as the mean projected distance between two clus- 
ters in the centre of the cluster is smaller than the merger radius 
(i.e. high filling factors) the merging process is very fast, because 
the clusters in the centre form a central merging object almost im- 
mediately and the subsequent merging happens preferably with this 
object which covers a fraction 77 of the surface area of the star form- 
ing region. Therefore the sub-clumps merge following a 'fast' ex- 
ponential decrease in their numbers: 

dr 



N(t) = iVo exp(-r/r). 



(5) 



77 is the fraction of the surface of the star-forming region, which is 
covered by the merger object and leads to a merging of the passing 
sub-clump. If time is measured in crosing times, i.e. the time the 
sub-clumps need to travel once through the star forming region, 
the merger rate can be derived by simple geometrical arg uments. 

~" d2002l). which 



We note that 77 replaces the e from iFellhauer et al 



we now reserve for the star formation efficiency. 
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Table 2. Values for &(a)No for the different choices of the parameters used in this study. The simulations are ordered in the same fashion as in Tabled] The 
first three columns show the parameters of the simulations, the fourth column the predicted values for 8(a)No and the fifth and sixth column show the values 
derived if we fit the theoretical formulae directly to the data of our simulations. The errors stated in those columns are the standard deviation of the fitting 
values and do not represent the much larger errors due to the low number statistics of the simulations. The final column gives the time needed to merge half of 
the sub-clumps f]/2i calculated using the fitted value from the slow regime. The time is given in real time, i.e. kyr. 
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d2002h find that t)y = 0.2 in their models. This value was found to 
be independent of the choice of a for a wide range of values. 

If the concentration of sub-clumps sinks below a threshold, 
i.e. when the mean projected distance between two clusters in the 
centre is larger than the merger radius, the 'slow' regime takes over. 
The sub-clumps still merge fast but rather with each other and build 
up the final merger object at later stages. Then the merging follows 



dN 
N(r) = 



—N 
No 



1 +5(a)N T 



(6) 



two clumps could have to merge after the interaction, which (mea- 
sured in units of the Plummer radius of the star for ming region) is 
a func tion of a. If we follow the same arguments as lFellhauer et al.l 
d2002l) . i.e. we start out fro m the energy cr iterion in the impulsive 
approximation developed bv lSpitzerl 111958): 
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and insert again t> tyP ~ %/2o" p i, r c = 3/2i? p i and take into 
account that the two clu sters will get gravitationally focussed 
dBinnev & Tremainefi 987) we arrive at the following expression: 
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which gives with 7 = 4.0 ± 0.5 the values for 8(a)No shown in 
Tab. [2] 7 is the mean of the projected radii of the outermost sub- 
clump in the realisations of the initial conditions, in units of R^{. 

The change between the very fast and not so fast regime takes 
place at a low a-value of < 0.02. 

Our merger criterion (i.e. when do we regard two sub-clumps 
as merged or a sub-clump with the merger object) is set up in the 
following way: if two clumps collide with each other and their cen- 
tres of density do not separate more than a certain radius r max 
for the rest of the simulation, then we regard the two clusters as 
merged. We favour a distance criterion to an energy criterion be- 
cause all clusters are bound to each other to begin with and in later 
stages we already have a merger object in the centre, which com- 
plicates the process. This merger object, which is somewhat larger 
than the initial clusters is also the reason to choose a rather large 
r mal . From hundreds of simulations we deduce that r max = R cu t, 
where R cu t is the cut-off radius of our single clumps, is a good 
choice for identifying two clumps as merged. It is small enough that 
the two single clumps do not separate from each other anymore and 
large enough to account for merger events with the central merger 
object as well. 



4 RESULTS 

4.1 Standard Model 

We begin by discussion the results for our standard model. In Fig.Q] 
we show how the number of s tar clusters dec r eases with time. 
According to the prediction of Fellhauer et al. ( 2002}) the num- 
ber should decrease exponentially with time (a > 0.02). The 
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Figure 1. Our standard model. Data points with error-bars show the mean 
values derived out of 5 simulations with different random seeds. The results 
of each of the simulations are shown as dotted lines. Dashed line (red) shows 
the analytical curve according to a 'fast' exponential decrease, solid (green) 
line shows the analytical curve according to the a-dependent 'slow' theory 
and dot-dashed (blue) line shows the fitted line according to an exponential 
decrease. 



data points can be nicely fit with an exponential decrease with 
rjat = 0.076 ± 0.003 with a reduced \ 2 of 1.00 (dot-dashed 
(blue) line in Fig. [T). This value is lower than the prediction of 

i c t = 0.2, i.e. the merging happens more slowly than pre- 
dicted. If we use instead the prediction of the 'slow' regime we 
obtain the solid (green) line in Fig.Q] This line seems to fit the data 
also. If we determine the fitted value of <5iVo for our data, we get 
(8N )fit = 0.10 ± 0.01 with a \ 2 of 2.78. This value is in very 
good agreement with the prediction of the a-dependent theory of 
(5A r o) P rcdict = 0.094. It seems that the prediction of the old the- 
ory works well, if we shift the transition value of a from the 'slow' 
to the 'fast' regime to a higher value. And if we check the velocity 
criterion formula (Eq. |4) we see (with our choice of No and e) that 
a ^ 0.07 ensures that the sub-clumps merge. 

Going back from our dimensionless time r (T^ — 260 kyr) 
to the real time t, we see that the majority of sub-clumps merge 
within 2-3 Myr, i.e. faster than any supernova explosion will oc- 
cur and blow out the remaining gas from the star-forming region. 
The model shows that by that time one already has a bound, nearly 
spherical merger object in the centre of the star forming region. 

We now have to find out why, and where, we have to switch 
from the very fast a-independent regime to the slower (but still 
fast) a-dependent theory. This is also interesting in the light of low 
mass star clusters with 10 3 to 10 4 stars, in which the gas removal 
happens even faster through stellar feedback of O stars. An example 
of such an object would be the Orion Nebula Cluster which has an 
age of 0. 5 to 1.5 Myr and has already blo wn out most of its natal 
material (Kroup a. Aarseth & Hur ley 2001). We therefore perform 
a suite of simulations with varying a, No and e. 



4.2 Parameter Study 

In Figs.[2]and[3]we show the results for the parameter survey. In the 
left panels of Fig.|2]we varied the value of a. It shows that as long 
as a ^ 0.2 the slow regime fits the data points well. Only if a be- 
comes very large (w 0.2) does the faster exponential decrease seem 
to take over. Finally, if the filling factor is as large as a = 0.5, most 
of the clumps overlap at the start of the simulation and the clumps 
do not separate but stick together and form a cluster from the start 
(small (blue) crosses in the top left panel). This general behaviour 
can be understood in terms of our velocity criterion, which tells us 
that as long as the filling factor is smaller than a « 0.07 merging 
between the clumps is possible. But by the time this threshold is 
reached, the filling factor is high enough for the clumps to form a 
central merger object from the start and the 'fast' regime takes over. 

If we keep a and e constant and only vary No the merging 
should happen according to the 'slow' regime with only a slight 
dependency on No. The right panels of Fig. [2] show this behaviour 
nicely, except for the bottom right panel, when iVo = 32. Here 
it seems that the clumps do not merge at all. Again this becomes 
clear with a look at Eq.[4] Inserting the constant values of e and a 
it results that No ^ 22 to allow the clumps to merge. This means 
by keeping all parameters constant and dividing the mass in stars 
into more and more clumps we reach the point where we leave the 
regime that two interacting clumps can absorb enough energy from 
the orbital motion that they finally merge. Merging in this regime 
is well suppressed and processes like destruction of the clumps due 
to the interactions take over. 

Finally we keep a and No constant and vary the star forma- 
tion efficiency e. The results are shown in Fig.f3]for two values of 
a. Using the fitted values for 8 No from Tab. [2] the data suggest a 
strong dependency of (Wo on e: 5 No oc e 2 . But a closer look re- 
veals that we rather have to deal with three different regimes. The 
first regime is governed by a strong background potential and al- 
most no clumps merging. Again the mass of a single clump is too 
low for the clump to absorb enough energy from the orbital motion 
into internal energy for merging to happen, i.e. the velocity crite- 
rion fails. If, in the a — 0.05 case (left panels), the star formation 
efficiency is higher than e > 0.25, the 'slow' a-dependent theory 
takes over and fits the data nicely. This again supports our velocity 
criterion which tells us that e ^ 0.24 is required for the clumps 
to be massive enough to allow merging. If the SFE is higher than 
approximately e = 0.5 then the 'fast' regime takes over, i.e. the 
central clumps are able to form a merger object immediately and 
clumps merge with the central object. 

In the right panels of Fig. [3] we see the same experiments for 
a higher value of a = 0.1. On a first glance one should expect 
that in this case the merging should happen faster, but the velocity 
criterion tells a different story. If we increase the filling factor by 
reducing the size of the star forming region (by keeping all other 
quantities constant) we raise the velocity dispersion of the clumps 
in the region and therefore make it more difficult for them to merge. 
Our simple velocity criterion (Eq. [4]l gives us a transition value of 
e ^ 0.47 while in our simulations we see that this transition has to 
be in the region of e sa 0.3. Also the 'fast' regime is reached only 
at a higher SFE- value of e ^ 0.7. 

As result of our parameter study we deduce tha t even t hough 
the simple geometrical theory of iFellhauer et al.l (2002) does 
not account for any background potential it is still applicable, 
with some lim i tation s. The velocity criterion of Eq. [3] used in 
Fellh auer et al] fc002h gives (translated into the values important 
in this study - Eq. [4]l a good hint for which part of the parame- 
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Figure 2. Left: Results of our simulations with varying values of ex. Solid (green) line shows the predicted decrease according to the 'slow' regime, while the 
dashed (red) curve gives the prediction of the 'fast' regime. The simulations follow the 'slow' regime except for very high values of a ^ 0.2. In the case of 
a = 0.5 (small (blue) crosses) in the top left panel all clumps overlap from the beginning and merge immediately (single simulation). The light (blue) crosses 
in lower right panel show a single simulation with a = 0.01, which also follow the slow regime. Right: Results of our simulations with varying Nq. The 
results show that the decrease in numbers is fairly independent of No and that the 'slow' regime (solid (green) lines; dashed (red) lines show the 'fast' regime) 
fits the data reasonably well (simulations with only 4 clumps suffer from 'very low number statistics') as long as the number of clumps are below a certain 
threshold which is reached at No = 32. See main text for explanation. In both sets of panels dotted lines are the results of the individual simulations. 
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Figure 3. Results of our simulations with different star formation efficiencies. Left: simulations with a = 0.05. Right: simulations with a = 0.1. The data 
points are the mean number of clusters after r crossing times out of all simulations with the same parameters. Solid (green) line is the analytical curve of 
the 'slow' regime, dashed (red) line is the analytical curve for the fast regime. Dotted lines represent each of the single simulations. One clearly sees that by 
varying e the results show three different regimes of the merging behaviour (see main text for explanation). 
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ter space the theory, especially the 'slow' regime, is valid. It still 
needs to be established when the transition between the 'slow' and 
the 'fast' regime sets in. From our parameter study, we can con- 
clude that a and e are anti-proportional to each other. However, the 
derivation of a general theory to explain the precise locations of the 
transition values of the parameters is postponed to a much wider 
study of the parameter space. 

To connect our theoretical study with actual observations we 
included in the last column of Tab.|2]the real time take for half of the 
sub-clumps to merge, which results already in a prominent merger 
object, taking the values of Tab. Q] This shows that for our set of 
parameters for the star forming cloud one can indeed find config- 
urations that form a merger object even before the stars manage to 
expell the gas via stellar winds and very long before the first super- 
nova is expected. More generally, the way to deduce this merging 
time for any given initial conditions is first to calculate 8 No: the 
'half-merger time' in crossing times of the star forming region is 
just the inverse of this number. Multiplying this number with the 
crossing time in years gives then the final answer. 

As a final remark one should keep the statistical nature of our 
results in mind. We use only 16 clumps to represent the distribution 
of clumps in the star forming region. In such random realisations 
it is very likely that even though we use a large filling factor, i.e. a 
high central concentration of clumps, the random numbers give us 
a realisation where most of the clumps are outside the central area 
and therefore act like a distribution with a much lower value of the 
filling factor. Also the opposite can be true. Furthermore by repeat- 
ing the experiment with the same set of parameters only 3-5 times 
we get closer to a mean representation of our input parameters but 
we also introduce extra noise to our study. In the figure-panels of 
this manuscript we show the behaviour of the single simulations as 
thin dotted lines to illustrate how differently random realisations of 
the same parameters can behave. 

4.3 Mass Segregation 

We mark the 50 innermost particles (< 1R P \) with the lowest en- 
ergy in each sub-clump and locate them in the final merger object. 
As expected they appear to end up in the centre of the merger ob- 
ject (see Fig. [4] top panel). To check if this is not merely a visual 
appearance due to the steep profile of the merger object, we derived 
the ratio of the surface density of the central particles to the surface 
density of all particles. As shown in the middle panel of Fig.|4]this 
ratio (£; here normalised to unity in the centre) declines as a shal- 
low power-law with index n = —0.56 ± 0.04. This shows that 
indeed mass-segregated stars in the centre of the sub-clumps end 
up again in the centre of the merger object. 

This process is a collective effect based on phase-space argu- 
ments only and is not due to two-body relaxation effects. Using 
a particle-mesh code for our simulations means that the effects of 
two-body relaxation are almost completely switched off. The signal 
is smeared out, due to the fact that we only choose particles accord- 
ing to their position and not simultaneously if they have a velocity 
which keeps them in the centre, which is the case with massive, seg- 
regated stars and disruptive effects of close encounters of unmerged 
sub-clumps with the merger-object. This means that the effect that 
central particles of sub-clumps end up in the centre of the merger 
object is not a process which requires two-body relaxation but 
rather a phase space argument. Still the results would be more pro- 
nounced if we would have used a direct N-body code, allowing for 
two-body relaxati on. This result nicely confirms and even s t rengfh - 
ens the results of McMillan, Vesperini & Porteg ies Zwartl J2007t) . 
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Figure 5. Lagrangian radii of one of our standard simulations after remov- 
ing the background potential (gas) instantaneously at the start of the simu- 
lation. One clearly sees that about 40 to 50 per cent of the cluster remains 
bound. 

who found that mass-segregated clumps form a mass-segregated 
star cluster. 



4.4 Infant Mortality 

Infant mortality is one of the major problems in forming star clus- 
ters: How can star clusters survive their gas-expulsion phase with- 
out getting completely destroyed? Theoretical work predicts that 
the star formation efficiency has to be e > 0.33 if the gas is lost 
instantaneously to retain a small bound core of a surviving clus- 
ter. There are several effects which can help to reduce this transi- 
tion value and therefore help star clusters to survive their formation 
process. 

First, if one takes into account internal evolution by simu- 
lating the cluster with an actual initial mass function with heavy 
stars sinking very fast to the centre, leading to a c ompact core, this 
value c an be reduced slightly (see experiments in lBoilv &~K roupa 
2003a. b). Furthermore it takes some time until the gas is removed 
from the cluster and the stars adapt to the change in potential. This 
helps especially very small star clusters (without massive stars) and 
very massive clusters (long gas-expulsion time). If star clusters do 
not form in isolation but in a star cluster complex the merging of 
the star clusters also helps to retain a massive star cluster in the 
form of the merger object. 

Here we add another mechanism which helps to overcome in- 
fant mortality. The time-scale of the merging in our standard model 
is very fast (even though it is the 'slow' regime). Within 2-3 Myr 
we see a merger object (star cluster) in the centre. Because the 
clumps have formed all over the star forming region and are now 
together in the centre, while the gas is decoupled from this merging 
process, the 'effective' star formation efficiency of the merged star 
cluster, i.e. the ratio between stars and total mass inside the merger 
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Figure 4. Left: Contour plot of the final merged cluster of one of our simulations. Over-plotted as (green) dots are the positions of the 50 initially lowest 
energy "marked" stars of each clump. Right: The ratio £ (crosses) between the surface-density profile of the marked stars and all stars of the merger object 
(normalised to 1) with the solid (red) line showing the power-law fit. The dashed (green) lines show the two profiles (also normalised to 1 in the centre) with 
the short dashed line showing all particles and the long dashed line the lowest (initially) energy particles only. 



object (not to be confused with t he effective SFE as usually found 
in the literature Goodwin (2008)), is higher than the overall effi- 
ciency of the parent cloud. The star cluster should therefore survive 
instantaneous gas-expulsion better. 

To show this we took the end-point of one of our merger simu- 
lations and removed the background potential instantaneously. The 
overall star formation efficiency of 0.32 should not allow the star 
cluster to survive, especially if modelled with a code which neglects 
internal evolution. But Fig.[5]shows that on the contrary we remain 
with a bound star cluster containing about 40 to 50 per cent of its 
initial mass. This can be explained by looking at the stars-to-gas 
ratio within the half-mass radius of our merger object. Here we see 
that by keeping the background (gas) potential constant the 'effec- 
tive' SFE_would have been of the order of 44 per cent (compare 
with lAdamsl ( feOOOh ). 

One can now argue that the gravitational pull of the newly 
formed star cluster would also collect some of the remaining gas. 
But we already start with a very simplistic realisation of the back- 
ground potential which is a strongly centrally concentrated Plum- 
mer model, allowing us to access all quantities analytically, being 
much more concentrated than a star forming region would be. So 
we indeed fulfill this criterion from the start. 



5 CONCLUSIONS 

In this study we tried to find out on which time-scales stars form- 
ing in clumps inside a star forming region merge and form an em- 
bedded star c luster. We have s hown that the analytical description 
developed by Fellha uer et alj J2002l) for the larger scale of merg- 
ing star clusters inside a cluster complex can be adopted for this 
problem also. 

By introducing a background potential, mimicking the re- 



maining gas, the basic theory need not be changed at all, although 
a shift of the critical values of the parameters which separate par- 
ticular merging behaviours is required. In the old study already a 
low filling factor of a 0.02 allowed the star clusters to merge 
very fast, decreasing their number exponentially. The presence of 
a background potential shifts this limit to rather high a values of 
about 0.2 due to the larger cluster-cluster velocities. This can be 
understood in the sense that the theory itself is geometrical and 
velocity independent, because all velocities were within a certain 
range which allows for this fast merging. So at each crossing time 
the same fraction of remaining clusters/clumps are able to merge 
with the central object. In the presence of a background potential 
the relative velocities of the clusters/clumps are higher, although 
the clusters have the same ability to absorb energy from the fly- 
bys. At any given time, this therefore reduces the number of clus- 
ters with velocities suitable for merging. 

Nevertheless, the old theory for the very low a regime (i.e. 
when the filling factor is so low that merging preferentially happens 
between two clusters/clumps at a time), provides a very good fit to 
the data. In this regime the number of clusters decreases according 
to 1/1 + t. For our standard model this translates into a merg- 
ing time for almost all of the sub-clumps which is much shorter 
than the time-scale until the first supernova explodes. Under these 
conditions a virialised star cluster forms well before gas-expulsion, 
driven by supernovae, occurs. For intermediate-mass clusters con- 
taining between about 10 3 and 10 4 stars gas expulsion occurs ear- 
lier and the requirement that sub-clustering disappears by the time 
the gas is expelled implies stringent constraints on the configura- 
tion of the sub-clusters. Nevertheless, such configurations are found 
in our parameter study. The new, hierarchically formed, embedded 
star cluster is then able to survive the gas-expulsion much more eas- 
ily than a monolithically formed cluster would have done, because 
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the 'effective' SFE at its location is much higher than the overall 
SFE of the region in which it formed. 

Finally this study has confirmed and st rengthened the findings 
of iMcMillan. Vesperini & Portegies Zwartl d2007» that mass- 
segregated clumps form mass-segregated clusters. Even though our 
method does not allow for mass-segregation at all, we are able to 
flag the lowest energy (i.e. the most bound) particles in the centre 
of each clump and locate them in the final merger object. We 
show that these particles are more centrally concentrated than the 
other particles. This shows that this process is governed by simple 
phase-space arguments and is not solely due to two-body relax- 
ation. As we have shown, it is not an effect of mass-segregation 
that massive, already mass-segregated (within their native clump) 
stars end up in the central region of the newly formed cluster. 
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